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Abstract 

The  integration  of  the  endothermic  reforming  reactions  and  the  exothermic  electrochemical  reactions  in  a  molten  carbonate  fuel  cell 
(MCFC)  is  investigated  with  a  mathematical  steady  state  anode  model.  The  dimensionless  model  consists  of  only  a  few  equations  based  on 
the  description  of  physical  phenomena,  and  is  thus  easily  solved  and  interpreted.  The  equations  and  their  derivation  are  presented  and  three 
different  applications  of  the  model  are  demonstrated.  It  proves  to  be  a  useful  tool  for  an  intuitive  understanding  of  the  interactions  of  reforming 
and  electrochemical  processes,  but  it  can  also  yield  quantitative  results  for  flow  sheet  design  of  MCFC  with  internal  reforming. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  storage  and  preparation  of  fuels  for  fuel  cells  is  cur¬ 
rently  under  investigation  by  numerous  research  and  develop¬ 
ment  groups  worldwide.  Its  impact  on  the  electric  efficiency 
and  profitability  of  any  type  of  fuel  cell  system  probably  is 
as  important  as  the  fuel  cell  design  itself.  Among  the  many 
proposed  and  investigated  variants  to  provide  hydrogen  for 
the  electrochemical  process,  the  steam  reforming  concept  has 
emerged  as  one  of  the  most  promising.  Fuels,  which  can  ei¬ 
ther  be  stored  easily  or  which  is  otherwise  available  at  the 
fuel  cell  site,  for  example  gasoline,  diesel  fuel  or  natural  gas, 
undergoes  a  reaction  which  produces  hydrogen  and  several 
by-products.  Often  times  this  reforming  product  has  to  be 
cleaned  up  before  being  fed  into  the  fuel  cell. 

The  key  issue  with  reforming  is  the  supply  of  heat  to  the 
endothermic  process  and  the  demand  for  high  temperatures 
in  order  to  obtain  high  reaction  rates  and  favourable  chemi¬ 
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cal  equilibria.  While  for  low  temperature  fuel  cells  this  heat 
has  to  be  provided  by  external  sources,  for  example  the  di¬ 
rect  oxidation  of  a  part  of  the  feed  gas,  high  temperature  fuel 
cells  like  the  SOFC  and  the  MCFC  directly  provide  the  heat 
required  for  the  reforming  process.  This  leads  to  the  inter¬ 
nal  reforming  concept,  in  which  the  endothermic  reforming 
process  is  integrated  into  the  fuel  cell,  where  the  exothermic 
electrochemical  reactions  occur.  The  interactions  of  both  pro¬ 
cesses  are  based  on  their  substantial  and  energetical  coupling 
and  their  effects  are  manifold.  The  use  of  mathematical  mod¬ 
els  is  probably  the  best  way  to  gain  deeper  understanding  of 
these  coherences. 

Some  authors  working  on  the  modelling  of  MCFC  include 
the  reforming  reactions  in  their  equations  in  different  ways. 
He  and  Chen  [1]  and  Yoshiba  et  al.  [2]  only  consider  the 
water-gas  shift  reaction  in  their  spacially  distributed  anode 
channel.  Due  to  its  high  rate  they  assume  it  to  be  in  thermo¬ 
dynamic  equilibrium.  Lukas  et  al.  [3]  and  Park  et  al.  [4]  also 
describe  the  water-gas  shift  reaction  as  being  in  equilibrium, 
but  they  additionally  include  the  steam  reforming  reaction 
of  methane  as  an  irreversible  reaction  with  a  finite  reaction 
rate  expression.  Especially,  Park  focuses  on  the  effect  of  an 
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Nomenclature 

Latin  symbols 

Arr j  Arrhenius  number  of  reaction  j 

Da;  Damkohler  number  of  reaction  j 

F  stoichiometric  factor 

kj  Arrhenius  term  of  reaction  j 

Kj  equilibrium  constant  of  reaction  j 

rj  reaction  rate  of  reaction  j 

Pel  electric  power 

R  anode  gas  recycle  ratio 

[/cell  cell  voltage 

Greek  symbols 

y  molar  flow  density  in  channel 

r  molar  flow 

r  corrected  molar  flow 

£  spacial  parameter 

iZ  temperature 

a air  ah'  number 

Vij  stoichiometric  coefficient  of  component  i  in  re¬ 

action  j 

v:  sum  of  stoichiometric  coefficients  of  reaction  j 

extent  of  reaction  j 
t-j' max  normalisation  factor  for 

<pE  electric  potential  in  electrolyte 

A (pjfi  standard  equilibrium  potential  difference  of  re¬ 

action  j 

Xi  molar  fraction  of  component  i 

Subscripts 

ref  reforming  reaction 

ox  oxidation  reaction 

red  reduction  reaction 

i  general  component  index 

feed  system  feed 

j  general  reaction  index 

max  see  above  §/,max 

0  standard  condition 

in  channel  inlet 

back  anode  gas  recycle 

out  channel  outlet 

ex  exhaust,  catalytic  burner  feed 

c  cathode 

air  burner  air 

Superscripts 

+  forward  reaction  direction 

—  backward  reaction  direction 

0  equilibrium 


indirect  reforming  unit  in  the  stack.  Although  all  of  these 
models  consider  the  influence  of  non-electrochemical  reac¬ 
tions  inside  the  fuel  cell,  none  of  them  is  suitable  to  easily  and 
intuitively  explain  their  impact  on  the  system.  This  is  mostly 
due  to  the  general  complexity  of  these  detailed  models. 

Other  authors  discuss  suitable  system  design  for  high  tem¬ 
perature  fuel  cells  with  internal  reforming.  Fellows  [5]  and 
Kortbeek  and  ottervanger  [6]  propose  complex  configura¬ 
tions,  including  gas  recycling,  side  feeding  and  cell  cascades. 
The  effects  of  these  measures  can  be  simulated  by  sophisti¬ 
cated  models,  but  these  do  not  offer  an  intuitive  understanding 
of  the  interaction  of  reforming  and  electrochemical  reactions 
nor  a  quick  and  simple  tool  for  conceptual  system  design. 

Therefore,  a  model  is  required  that  describes  the  basics  of  a 
high  temperature  fuel  cell  in  a  physically  meaningful  way  and 
fulfils  these  requirement.  In  addition,  an  easily  interpretable 
diagram  is  useful  for  the  illustration  of  the  numerical  results. 
The  objective  of  this  paper  is  to  develop  such  a  tool  on  the 
example  of  a  methane-fed  MCFC,  and  to  demonstrate  its 
usefulness  for  several  applications. 


2.  Technical  background 

The  working  principle  of  the  Molten  Carbonate  Fuel  Cell 
is  illustrated  in  Fig.  1.  The  anode  is  fed  with  a  preheated 
mixture  of  desulphurised  natural  gas  and  steam  at  a  steam  to 
carbon  ratio  of  about  S/C  =  2.5.  This  feed  gas  is  converted 
via  steam  reforming  into  a  hydrogen  rich  gas  mixture  at  the 
reformer  catalyst,  which  is  placed  inside  the  anode  channel. 
Carbon  monoxide  is  a  major  by-product  of  this  reforming 
process.  Simultaneously,  the  water  gas  shift  reaction  trans¬ 
forms  carbon  monoxide  into  carbon  dioxide  and  additional 
hydrogen: 

CH4  +  H2O  CO  +  3H2  (1) 

CO  +  H20  -o-  CO2  +  H2  (2) 

The  equilibrium  limitations  of  these  two  reforming  reac¬ 
tions  are  overcome  by  continuous  removal  of  hydrogen  and 
carbon  monoxide,  which  is  directly  oxidised  electrochemi- 
cally  at  the  anodic  electrode.  There,  these  substances  react 
with  carbonate  ions  from  the  electrolyte  producing  carbon 


CH4+H70<=>C0  +  3H,  H,0 
CO  +H20<=>C0+Hi  CO, 

H2  +co32  <=>h20+  co2+  2e 
CO+CO/  »  ‘  2CO+2e 


02 

N„ 


Anode  electrode 
Electrolyte 
Cathode  electrode 

Exhaust  gas*" 


T~600°C 


Catalytic 
Combustion 
Chamber 

T~ 

fo'<-co2 

Cathode  gas  recycle 


Fig.  1 .  Working  principle  of  the  MCFC  with  direct  internal  reforming. 
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dioxide,  water  and  electrons  according  to  the  following  stoi¬ 
chiometric  relations: 

H2  +  C032“  H20  +  C02  +  2e“  (3) 

CO  +  C032-  2C02  +  2e“  (4) 


In  addition  to  the  mass  coupling,  the  reforming  process, 
which  is  overall  endothermic,  is  energetically  coupled  to  the 
exothermic  oxidation  reactions,  thus  the  heat  required  for  the 
generation  of  hydrogen  is  provided. 

The  anode  exhaust  gas  is  mixed  with  air  and  the  unoxi¬ 
dised  components  are  fully  oxidised  in  a  catalytic  combustion 
chamber.  Because  air  is  fed  in  excess  the  exhaust  gas  from 
the  burner  still  contains  a  significant  amount  of  oxygen.  This 
gas  is  then  fed  to  the  cathode  channel  where  the  electrochem¬ 
ical  reduction  takes  place.  New  carbonate  ions  are  produced 
from  carbon  dioxide  and  oxygen  according  to  the  backward 
direction  of  the  following  cathode  reaction: 

C032“  \02  +  C02  +  2e"  (5) 

The  carbonate  ions  are  transported  towards  the  anode  elec¬ 
trode  through  the  electrolyte,  which  is  a  eutectic  carbonate 
melt.  Often  times,  a  part  of  the  cathodic  exhaust  gas  is  recy¬ 
cled  back  towards  the  catalytic  combustion  chamber  and  the 
remaining  gas  leaves  the  fuel  cell.  Due  to  its  high  temperature 
of  about  600  °C,  the  exhaust  gas  can  be  used  for  subsequent 
steam  generation,  for  additional  production  of  electric  energy 
via  a  micro  turbine  or  for  other  purposes. 

Besides  the  possibility  to  form  hydrogen  inside  the  an¬ 
ode  channel,  it  can  also  be  produced  otherwise.  Fig.  2  shows 
the  three  basic  reforming  concepts  for  high  temperature  fuel 
cells.  The  external  reforming  (ER)  happens  in  a  separate  ves¬ 
sel  outside  the  cell  stack.  Any  arbitrary  temperature  can  be 
applied,  but  high  temperatures  require  external  heating.  For 
temperatures  well  below  the  cell  temperature,  the  cell  exhaust 
gas  can  be  used  to  supply  the  required  reaction  heat  to  the 
ER.  The  main  advantage  of  the  ER  is  that  its  design  and  op¬ 
erating  conditions  are  largely  independent  from  those  of  the 
fuel  cell. 

The  indirect  internal  reformer  (HR)  is  situated  in  the  cell 
stack  in  special  reforming  channels,  where  only  reforming 
reaction  takes  place.  This  concept  features  energetic  coupling 
with  the  exothermic  oxidation  process.  The  main  advantage  is 
that  no  external  heat  exchanger  is  required,  as  the  separator 
plate  between  HR  and  anode  channel  adopts  this  function. 


T=T 

ex 


Air 


ER:  External  Reforming 
HR:  Indirect  Internal  Reforming 
DIR:  Direct  Internal  Reforming 


Fig.  2.  Reforming  concepts  for  high  temperature  fuel  cells. 


The  HR  can  be  seen  as  an  external  reformer  operating  at  fuel 
cell  temperature. 

Finally,  the  substantial  and  energetic  coupling  of  the  re¬ 
forming  and  electrochemical  processes  leads  to  the  direct  in¬ 
ternal  reforming  (DIR)  concept.  Not  only  is  heat  exchanged 
between  both,  but  also  the  products  of  the  reforming  process 
are  continuously  removed  by  the  oxidation  reactions.  This  is 
a  highly  integrated  concept,  which  offers  the  chance  for  high 
system  efficiency,  but  it  requires  sophisticated  system  design 
in  order  to  work  properly. 

3.  Modelling 

With  the  numerous  possible  combinations  of  different  re¬ 
forming  concepts  and  due  to  the  various  flow  sheet  options 
like  gas  recycling,  cell  cascading  and  side  feeding,  a  tool  is 
required  for  the  analysis  of  these  configurations.  This  tool 
should  be 

•  physically  meaningful  in  order  to  allow  for  an  intuitive 
understanding  of  its  results, 

•  flexible,  so  that  one  can  also  analyse  additional  aspects 
besides  those  mentioned  here, 

•  fast  to  solve,  as  the  task  here  is  not  a  detailed  and  exact  pre¬ 
diction  of  states,  but  should  rather  yield  principal  insight 
and  show  up  tendencies, 

•  easily  illustrated  and  interpretable. 

In  this  section,  a  mathematical  model  called  the  steady 
state  anode  model  is  introduced  which  fulfils  these  require¬ 
ments.  First,  the  complete  model  equations  are  given  and  their 
derivation  is  briefly  indicated.  In  the  second  part  an  accom¬ 
panying  diagram  to  this  model  is  introduced,  in  which  the 
simulation  results  can  be  displayed. 

3.1.  Model  derivation 

The  steady  state  anode  model  is  derived  from  a  fully 
dynamic,  spacially  two-dimensional  description  of  a  single 
cross-flow  MCFC  [7,8].  As  thar  model  is  completely  formu¬ 
lated  in  terms  of  dimensionless  variables,  the  steady  state 
anode  model  is  dimensionless  as  well.  It  is  derived  under  the 
following  assumptions: 

•  The  system  is  at  steady  state. 

•  Assume  a  sufficiently  high  cathode  gas  recycle  ratio,  so 
that  the  composition  and  temperature  of  the  cathode  gas 
are  approximately  constant  along  the  channel. 

•  As  a  consequence  of  the  constant  cathode  gas  temperature, 
the  cell  temperature  is  also  constant  along  the  spacial  co¬ 
ordinate.  This  temperature  is  not  calculated  by  an  energy 
balance,  but  it  is  given  as  a  parameter.  Nevertheless,  tem¬ 
perature  effects  on  the  reaction  rates  and  the  equilibrium 
constants  are  considered. 

•  The  cell  is  operated  under  ambient  pressure  conditions. 
Pressure  drop  is  neglected,  as  we  mainly  focus  on  the  sub¬ 
stantial  interplay  of  reforming  and  oxidation. 
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•  The  oxidation  reaction  of  carbon  monoxide  (Eq.  (3))  in 
MCFC  is  known  to  be  of  minor  importance  compared  to 
the  oxidation  of  hydrogen  and  is  thus  neglected. 

•  Due  to  the  very  fast  water-gas  shift  reaction  and  because 
the  equilibrium  of  this  reaction  at  high  temperatures  is  on 
the  right  hand  side  of  Eq.  (2),  the  amount  of  carbon  monox¬ 
ide  in  the  anode  gas  is  negligible.  Thus,  both  reforming  re¬ 
actions  are  condensed  into  one  single  reforming  reaction: 

CH4  +  2H20  C02  +  4H2  (6) 

•  Mass  transport  in  the  electrode  pores  is  neglected. 

•  Instead  of  galvanostatic  cell  operation,  potentiostatic  op¬ 
erating  mode  is  applied.  Thus,  the  cell  voltage  is  given  as 
an  operating  parameter. 

•  The  ion  conduction  resistance  in  the  electrolyte  perpen¬ 
dicular  to  the  cell  plain  is  negligible.  Nevertheless,  no  ion 
conduction  occurs  along  the  cell  plain.  Thus  the  electrolyte 
has  identical  electric  potentials  at  the  anode  and  cathode 
electrode,  but  these  are  locally  distributed. 

•  The  possibility  of  anode  gas  recycling  shall  be  considered. 


So 


r  k 

back  ^back 


0  ' '  fin T 


u 


’  ^  Hou. 

— 

Catalytic 

Combustion 

L 

Chamber 

Fig.  3.  Flow  scheme  and  assigned  variables  of  the  steady  state  anode  model. 


with  this,  the  molar  fraction  of  any  component  can  be  de¬ 
scribed  as  a  function  of  the  extents  of  reaction: 

/'.feed  +  J2j=ref,oxviJ  x  §/£/,max 

Xi  =  — .  ^ - 1 - rr -  (7) 

1  +  2—, y=ref,oxk/  x  SjSj, max 


where 


XCH4,feed 


1 

1  +  S/C 


The  resulting  equations,  i.e.  mainly  the  partial  and  the  to¬ 
tal  mass  balances  in  the  anode  channels,  are  given  in  dimen¬ 
sionless  form  in  Appendix  A.  The  five  ordinary  differential 
equations  (ODEs)  describe  five  unknowns  in  one  single  an¬ 
ode  channel,  namely  the  molar  fractions  of  CH4,  H?0,  H2 
and  CO?  as  well  as  the  molar  flow  density,  y.  The  cathode 
gas  composition  is  calculated  independently  from  the  ODEs 
according  to  a  full  combustion  calculation  of  the  feed  stream 
with  air  in  several  algebraic  equations  (see  Appendix  A). 

To  reduce  the  number  of  ODEs  even  further,  a  physically 
motivated  transformation  is  applied.  The  idea  is  that  for  a 
given  composition  of  the  feed  gas  (CH4/H20-mixture,  char¬ 
acterised  by  the  steam/carbon-ratio)  the  composition  of  the 
gas  at  any  point  in  the  anode  channel  or  in  any  of  the  re¬ 
forming  units  is  described  by  only  two  states:  the  extent  of 
the  reforming  reaction,  assigned  £ref,  and  the  extent  of  the 
oxidation  reaction,  assigned  £ox.  These  variables  are  made 
dimensionless  and  normalised  to  unity,  so  they  can  be  inter¬ 
preted  as  follows: 

|ref  =  0  gas  is  not  reformed  . . . 

|ref  =  1  gas  is  completely  reformed 


|ox  =  0  gas  is  not  oxidised  . . . 

|ox  =  1  gas  is  completely  oxidised 


XHaO.feed 


S/C 

1+S/C 


/ref 


,max  — 


1 

1  +  S/C 


4 

lox.max  -  j  +S/C 

In  addition  to  this,  the  possibility  of  anode  gas  recycling 
shall  be  considered.  Fig.  3  shows  the  assignments  of  all  molar 
gas  flows  inside  and  outside  the  anode  channel.  At  steady 
state,  the  flows  differ  due  to  two  reasons:  firstly  the  addition 
of  the  recycle  flow  to  the  anode  inlet,  and  secondly  due  to  the 
change  of  total  mole  numbers  with  reactions. 

To  eliminate  the  influence  of  the  second  mentioned,  a  cor¬ 
rected  molar  flow  is  introduced,  marked  with  a  tilde  sign.  In 
general,  it  holds  that 


r  —  C  x  I  1  +  ^2  vj  x  Hj%j, max  I  (8) 

\  ;=ref,ox  / 

with  this  definition,  the  mixing  rules  can  be  set  up.  They 
follow  from  the  total  and  the  partial  molar  balance  at  the 
mixing  point  before  the  channel  inlet  (Fig.  3)  together  with 
Eqs.  (7)  and  (8).  The  flows  are  added,  while  the  mixing  gas 
composition  is  a  linear  combination  of  the  gas  compositions 
of  both  mixed  flows  according  to  the  law  of  the  lever: 


note  that  the  extent  of  oxidation  reaction  is  identical  to  the 
well  known  fuel  utilisation,  relating  the  rate  of  electrochem¬ 
ical  fuel  consumption  to  the  fuel  feed  rate.  Usually,  a  high 
fuel  utilisation  is  desirable. 

|ox+fuel  utilisation 


f  in  —  /()  +  fback 


(9) 


Hj,  in  — 


Cq$  ;,0  +  /b  ack/  /.back 

To  +  Tback 


(10) 


The  extents  of  reaction  inside  the  channel  are  described  by 
ordinary  differential  equations  resulting  from  the  continuous 
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total  and  partial  molar  balances.  Their  derivation  is  given  in  The  second  equation  sums  the  potential  differences  at  both 
Appendix  A.  The  result  reads:  electrodes,  which  results  in  the  given  cell  voltage: 


=  D'dJri  x  1 
fin  Hj,  max 

with 


(ID 


j  =  ref,  ox 

The  boundary  condition  is  given  at  the  inlet 


Hjtt  =  0)  =  SjM  (12) 

The  Damkohler  numbers  for  the  reforming  and  the  oxi¬ 
dation  reaction  cannot  only  be  used  to  adjust  finite  reaction 
rate  constants  of  both  reactions  in  the  simulation,  but  each 
of  these  reactions  can  be  eliminated  by  setting  the  specific 
Damkohler  number  to  zero.  Thus  it  is  possible  to  simulate  a 
pure  reforming  channel  (Daox  =  0)  or  an  anode  channel  with¬ 
out  direct  internal  reforming  (Dareci  =  0). 

In  addition  to  the  two  ODEs,  the  reaction  kinetics  are  re¬ 
quired.  The  kinetic  expressions  are  given  here  in  short  form. 
For  a  full  version  of  these  equations  we  refer  the  reader  to 
Appendix  A.  The  reforming  reaction  rate  depends  on  the  gas 
composition  and  the  given  temperature  parameter,  the  equi¬ 
librium  constant  is  only  temperature  dependent: 


A-ef©  ft)  =  *refW  X  r+f©  -  rref®  (13) 

I  I 

The  oxidation  rate  depends  not  only  on  the  gas  composi¬ 
tion  and  the  temperature  parameter,  but  also  on  the  electric 
potential  difference  at  the  charged  double  layer  of  the  anode 
electrode,  Ar/>a.  On  the  other  hand,  the  reduction  reaction 
rate  depends  on  the  electric  potential  difference  at  the  cath¬ 
ode  electrode,  A 4>a.  Applying  a  Butler- Volmer  type  reaction 
kinetics,  we  receive: 


A>x©  ft  A</>a)  —  km(i).  Af/;a)  x 


x©~ 


x© 


^eq,ox(ft  ^</>a) 


(14) 


^red(£?  ft  A(/>c)  —  &ox(ft  A0C) 


X 


rred©  1 

A cq,  red  (ft-  A  </>c)  J 


(15) 


With  the  Eqs.  (14)  and  (15)  containing  four  unknowns, 
namely  rox,  rred,  A</;a  and  A <pc,  two  more  equations  are 
needed  to  solve  the  equations.  One  of  them  is  the  steady 
state  charge  balance  stating  that  the  current  density  produced 
at  the  anode  electrode  must  be  equal  to  the  current  density 
produced  at  the  cathode  electrode  at  any  given  point  in  the 
channel: 


Daox  x  rox  —  Dared  x  (  r red) 


(16) 


A0a  +  (-A0c)=  t/cell  (17) 

Eqs.  ( 1 4)— ( 17)  contain  four  unknowns:  the  anodic  and  ca¬ 
thodic  reaction  rates,  rox  and  rox,  respectively,  and  the  electric 
potentials  at  both  electrodes,  A <pa  and  A (pc.  They  have  to  be 
solved  for  the  oxidation  reaction  rate,  rox,  which  is  required 
inEq.  (11). 

As  a  measure  for  the  system  efficiency,  the  electric  power 
output  of  a  cell  is  calculated  according  to 

(18) 

F 

where 

Afox  =  fox©  =  fend)  -  fox©  =  0)  (19) 

With  this,  the  complete  model  is  described.  Outside  any  chan¬ 
nels,  Eqs.  (9)  and  (10)  describe  mixing  rules.  These  might  be 
necessary  to  calculate  the  boundary  condition  (Eq.  (12))  for 
Eq.  (11),  which  together  with  the  reaction  rate  expressions  in 
Eqs.  ( 13)— ( 17)  describe  the  state  inside  such  channels. 

3.2.  The  conversion  diagram 

With  only  two  locally  distributed  states  in  the  model,  it  is 
possible  to  illustrate  any  of  its  results  in  a  phase  diagram.  In 
this  kind  of  plot,  the  independent  spacial  variable  is  elimi¬ 
nated  by  plotting  one  state  versus  the  other.  Fig.  4  shows  the 
extent  of  the  reforming  reaction  on  the  horizontal  axis  and 
the  extent  of  oxidation  reaction  on  the  vertical  axis.  Because 
the  molar  fraction  of  any  component  cannot  become  negative, 
certain  limitations  for  physically  reasonable  states  in  the  con¬ 
version  diagram  follow.  The  extents  of  reaction  may  only  be 
situated  within  the  triangle  in  Fig.  4. 

Also  indicated  in  this  figure  is  the  fact  that  for  a  given 
temperature  and  cell  voltage,  the  reaction  rates  at  every  point 
within  the  allowable  region  are  defined.  A  positive  reforming 


Conversion  diagram 


Fig.  4.  The  conversion  diagram,  stoichioinetrically  attainable  region  and 
reaction  rate  vector. 
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Conversion  diagram  Conversion  diagram 


Fig.  5.  Reforming  equilibrium  lines  at  various  temperatures  and 
oxidation  equilibrium  lines  at  constant  temperature  and  var¬ 
ious  cell  voltage.  The  dimensionless  numbers  correspond  to 
the  following  values:  0  =  {2.5,  3  and  4}  =  {470,  620  and  920  °C} 
Uct\\  =  {31,  34  and  38}  =  {796,  873  and  976 mV). 

reaction  rate  means  that  the  gas  state  moves  to  the  right,  while 
a  positive  oxidation  rate  causes  a  move  upwards.  Together, 
both  reaction  rates  form  a  vector  indicating  in  which  direction 
the  state  moves  along  the  spacial  coordinate. 

All  reactions  are  considered  to  be  reversible,  so  the  gas 
composition  may  attain  equilibrium  for  one  or  both  reac¬ 
tions  in  the  anode  channel.  In  Fig.  5,  several  equilibrium 
lines  are  displayed  for  different  temperatures.  On  these  lines, 
the  reaction  rate  of  the  corresponding  reaction  equals  zero. 
The  reforming  equilibrium  lines  start  from  the  right  bot¬ 
tom  corner  and  end  in  the  right  upper  corner  of  the  ad¬ 
missible  region.  How  far  they  reach  to  the  left  only  de¬ 
pends  on  the  given  temperature.  In  general,  low  tempera¬ 
tures  come  with  a  low  equilibrium  conversion  of  the  re¬ 
forming  reaction,  thus  the  equilibrium  curve  reaches  far  to 
the  left.  Reforming  equilibrium  lines  can  be  calculated  from 
Eq.  (13)  by  setting  the  rate  to  zero  and  then  solving  for  a 
function  §ref(§ox),  at  which  the  reforming  reaction  rate  is 
zero. 

The  oxidation  equilibrium  line  always  starts  at  the  dia¬ 
grams  origin  point  and  end  somewhere  on  the  right  bound¬ 
ary  line.  Their  course  depends  on  the  given  parameters  of 
temperature  and  cell  voltage.  The  three  exemplary  lines  in 
Fig.  5  are  at  constant  temperature,  but  varying  cell  voltage. 
On  such  a  line,  the  oxidation  process  comes  to  a  halt.  The 
oxidation  equilibrium  line  is  calculated  from  Eqs.  (14),  (15), 
and  (17),  where  both  electrochemical  reaction  rates  are  set 
to  zero.  The  three  resulting  equations  then  contain  two  inde¬ 
pendent  variables,  which  are  both  extents  of  conversion,  and 
two  unknown  potential  differences.  This  allows  to  solve  for 
a  function  £0x(£ref)>  at  which  electrochemical  equilibrium  is 
given. 

Due  to  their  characteristic  course,  there  is  always  exactly 
one  intersection  point  of  both  curves  for  given  temperature 
and  voltage.  At  this  point,  both  reaction  rates  become  zero. 


Fig.  6.  Two  exemplary  trajectories  for  identical  temperature  and  cell  voltage, 
but  with  different  inlet  conditions.  Both  end  up  in  the  same  attractor  point, 
which  is  the  intersecting  point  of  the  reforming  and  the  oxidation  equilibrium 
lines. 

Thus  this  is  a  stationary  point.  Because  the  reactions  always 
run  towards  their  equilibrium,  this  stationary  point  is  an  at¬ 
tractor. 

Two  exemplary  trajectories  are  shown  in  Fig.  6.  The 
boundary  condition  of  the  first  simulation  is  that  both  ex¬ 
tents  of  reaction  are  zero,  meaning  that  the  gas  is  neither 
reformed  or  oxidised  the  channel  entrance.  Upon  entering, 
the  reforming  reaction  starts  and  as  soon  as  there  is  a  signif¬ 
icant  amount  of  hydrogen  available,  the  oxidation  reaction 
rate  also  increases.  Thus,  the  trajectory  moves  from  the  ori¬ 
gin  to  the  right  and  then  starts  to  move  upward  as  well.  For 
an  infinite  channel  length,  the  trajectory  ends  in  the  attractor 
point.  The  same  holds  true  for  a  different  inlet  condition.  As¬ 
sume  a  completely  reformed  gas  entering  the  anode  channel. 
It  ends  up  in  the  same  attractor  point  as  does  the  other  trajec¬ 
tory,  but  on  its  way,  the  reforming  reaction  runs  in  backward 
direction  and  the  reforming  equilibrium  line  is  crossed  in  a 
vertical  manner. 


4.  Applications 

The  steady  state  anode  model  can  be  applied  for  numerous 
studies.  In  the  following,  three  of  them  are  demonstrated. 
They  are 

•  a  comparison  of  different  combinations  of  reforming  con¬ 
cepts, 

•  the  simulation  and  numerical  evaluation  of  the  benefit  of 
a  fuel  cell  cascade, 

•  the  simulation  of  an  anode  exhaust  gas  recycle. 

4.1.  Comparison  of  reforming  concepts 

As  indicated  in  Fig.  2,  three  different  concepts  are  known 
for  reforming  in  high  temperature  fuel  cell  systems.  The 
steady  state  anode  model  allows  comparing  the  various  com- 
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Fig.  7.  Trajectories  for  various  reforming  configurations  without  DIR. 


binations  of  these  in  an  illustrative  way.  First,  consider  a  sys¬ 
tem  without  reforming  catalyst  inside  the  anode  channel,  that 
means  without  DIR.  Three  alternatives  for  the  fuel  gas  treat¬ 
ment  are  discussed: 

•  A  low  temperature  ER,  which  can  easily  be  heated  by  the 
hot  exhaust  gases  of  the  cell. 

•  An  HR,  where  the  reforming  temperature  is  identical  to  the 
cell  temperature. 

•  A  high  temperature  ER,  which  requires  additional  external 
heating. 

In  this  application,  all  three  reforming  units  are  considered 
to  be  infinitely  long,  thus  equilibrium  is  reached  in  each  unit. 
The  exhaust  gas  of  the  reforming  unit  is  fed  into  the  anode 
channel,  where  only  oxidation  reaction  takes  place. 

The  resulting  trajectories  are  plotted  in  Fig.  7.  Low  temper¬ 
ature  ER  and  HR  produce  anode  fuel  gas,  which  is  reformed 
to  an  extent  of  30  and  65%,  respectively.  The  following  fuel 
cell  can  only  reach  fuel  utilisations  of  25  and  60%,  which 
is  unsatisfying,  because  large  parts  of  the  fuel  are  not  used. 
Only  the  high  temperature  ER  offers  the  possibility  to  reach 
high  fuel  utilisation  of  about  90%,  but  its  clear  drawback  is 
the  need  for  extra  heating,  which  surely  would  decrease  the 
overall  efficiency.  Thus  none  of  these  alternatives  is  satisfy¬ 
ing. 

The  situation  changes  significantly  when  DIR  is  applied. 
Fig.  8  shows  that  independent  of  the  pre-treatment  of  the 
gas,  high  fuel  utilisations  can  be  reached,  as  all  trajectories 
end  in  one  single  point  close  to  the  right  upper  corner.  Here, 
the  high  temperature  ER  has  one  additional  drawback:  As 
the  cell  temperature  is  lower  than  the  reforming  temperature, 
the  reforming  reaction  runs  in  backward  direction  in  the  first 
part,  spoiling  the  high  degree  of  reforming  conversion.  The 
same  attractor  point  can  be  reached  by  a  system  without  pre¬ 
treatment,  where  the  fresh  gas  is  directly  fed  into  the  anode 
channel. 

This  consideration  clearly  shows  the  advantage  of  the  DIR 
concept.  Independent  of  the  pre-treatment  of  the  feed  gas, 
the  application  of  DIR  always  leads  to  high  degrees  of  fuel 


Fig.  8.  Trajectories  for  various  reforming  configurations  with  DIR.  High 
degrees  of  fuel  utilisation  are  reachable.  Compare  to  Fig.  7. 

utilisation.  The  steady  state  anode  model  together  with  its 
conversion  diagram  serves  as  an  illustrative  tool  here. 

4.2.  Anode  cascade 

Various  authors  consider  appropriate  flow  sheets  for  fuel 
cell  systems.  Fellows  [5]  proposes  to  employ  a  cascade  of 
fuel  cells,  operating  at  different  cell  voltages,  to  increase  the 
overall  electric  power  output.  The  proposal  is  that  the  exhaust 
gas  of  the  first  cells  anode  channel  is  fed  into  a  second  cell 
and  so  on.  The  cascade  discussed  here  consists  of  two  cells 
(Fig.  9). 

Concerning  electric  connections,  two  different  configura¬ 
tions  are  possible  (Fig.  10).  If  each  cell  is  operated  not  only 
at  its  own  cell  voltage,  but  it  also  has  an  independent  cell  cur¬ 
rent,  they  are  connected  as  parallel  current  sources  and  both 
require  a  separate  electrical  converter  (Fig.  10a).  The  alter¬ 
native  is  that  although  their  cell  voltage  differs,  the  same  cell 
current  runs  through  both  cells  (Fig.  10b).  This  series  con- 


Fig.  9.  Flow  sheet  for  fuel  cell  cascade. 


o- 


J=Mi 


(b) 


Fig.  10.  Electric  configurations  of  fuel  cell  cascade:  (a)  parallel  configura¬ 
tion  with  independent  cell  voltage  and  cell  current;  (b)  series  configuration 
with  independent  cell  voltage  and  equal  cell  current. 
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Fig.  1 1 .  Trajectories  for  fuel  cell  cascade  in  electric  parallel  configuration  Fig.  12.  Trajectories  for  fuel  cell  cascade  in  electric  series  configuration  (for 
(for  quantitative  results  see  Table  1).  quantitative  results  see  Table  1). 


nection  requires  only  one  electrical  converter,  thus  it  needs 
less  hardware  effort.  Both  versions  are  discussed  here. 

The  two  cascade  configurations  are  compared  to  a  single 
cell,  which  has  the  same  size  as  the  sum  of  the  cell  sizes  in 
each  cascade.  Each  system  is  optimised  to  yield  an  optimum 
of  electric  power  at  a  constant  feed  rate.  The  optimisation 
variables  are  the  cell  voltage  for  the  single  cell  system,  and 
both  cell  voltages  plus  the  size  of  the  first  cell  for  both  cascade 
configurations. 

Fig.  11  shows  the  resulting  trajectories  for  the  parallel 
configuration  in  comparison  to  the  single  cell.  The  electric 
power  output  is  illustrated  by  the  rectangles  on  the  left  of 
the  vertical  axis.  Their  area  is  the  product  of  the  cell  volt¬ 
age  and  the  advance  in  oxidation  conversion  over  each  cell. 
Table  1  shows  that  the  electric  power  output  of  the  cascade 
configuration  exceeds  that  of  the  single  cell  by  4.6%. 

The  trajectory  for  the  series  connection  is  shown  in  Fig.  12. 
The  differences  of  oxidation  conversion  over  both  cells  are 
equal,  so  the  additional  condition  of  equal  cell  currents  in  both 
cells  is  fulfilled.  This  condition  leads  to  a  smaller  increase 
of  electric  power,  nevertheless  an  additional  benefit  of  3.7% 
due  to  the  cascade  configuration  at  constant  feed  rate  can  be 
considered  relevant. 

According  to  these  calculations,  the  option  of  applying 
fuel  cell  cascades  to  increase  the  electric  power  output  and 
thereby  increase  the  system  efficiency  is  an  interesting  option. 
From  an  economical  point  of  view,  the  additional  hardware 


Table  1 

Optimisation  results  of  single  cell  and  both  two-cell  cascades 


Length  f 

Voltage 

f^cell 

Extent  of  oxidation 
reaction  A£0x 

El.  power  Pei 

Single 

1.0 

28.83 

0.947 

62.40 

Parallel  cascade 

1.  Cell 

0.674 

30.81 

0.683 

48.10 

2.  Cell 

0.326 

27.00 

0.278 

17.17 

Total 

1.0 

- 

0.961 

65.27  (+4.6%) 

Series  cascade 

1.  Cell 

0.508 

31.84 

0.480 

34.93 

2.  Cell 

0.492 

27.20 

0.480 

29.84 

Total 

1.0 

- 

0.960 

64.70  (+3.7%) 

effort  necessary  to  build  such  a  system  stands  against  this 
configuration,  at  least  for  small  systems.  For  large  systems, 
where  several  stacks  are  combined  in  a  plant,  this  might  be 
profitable.  The  benefit  of  the  steady  state  anode  model  in  this 
case  goes  beyond  illustration  and  intuitive  understanding,  but 
it  can  also  serve  here  as  a  tool  for  a  first  quantitative  system 
design  and  it  allows  conducting  first  economical  calculations. 


4.3.  Anode  exhaust  gas  recycling 


In  literature  one  can  also  find  the  proposal  to  recycle  a  part 
of  the  anode  exhaust  gas  back  towards  the  inlet  [5,6]  (Fig.  3). 
Especially  for  pure  DIR  systems  this  has  certain  advantages: 
if  fresh  feed  gas  is  fed  into  the  anode  channel,  no  hydrogen 
is  present  at  the  inlet  and  thus  the  electrode  is  unused  in  this 
region.  This  can  be  amended  by  recycling  part  of  the  exhaust 
gas,  which  still  contains  some  hydrogen,  back  to  the  inlet. 
Define  the  recycle  flow  via  a  recycle  ratio  according  to: 

R  =  5^-  (20) 

Tback 

Then  the  mixing  rules  (Eqs.  (9)  and  (10))  simplify  to 


^in  =  (1  +  R)  Peed 


(21) 


Hj,  in  — 


Hj,  0  %  back 

rr ~r 


(22) 


Fig.  13  shows  the  trajectory  of  a  single  cell  without  ( R  =  0) 
and  with  anode  gas  recycling  (R  =  0.3).  Both  are  operated 
at  individual  optimal  cell  voltage,  so  that  electric  power  is 
maximised.  While  the  inlet  condition  of  the  system  without 
recycle  is  at  the  origin,  the  inlet  condition  of  the  recycle  sys¬ 
tem  is  calculated  via  the  law  of  the  lever  (Eq.  (19))  from  the 
conditions  of  the  feed  gas  and  the  anode  exhaust  gas. 

The  electric  power  output  of  the  system  with  gas  recycle 
is  always  lower  than  what  one  gets  from  the  system  with¬ 
out  recycling.  This  can  be  explained  in  two  ways,  both  of 
them  have  the  same  basis.  The  first  explanation  is  that  back 
mixing  always  worsens  the  reactor  performance  in  systems 
with  a  positive  order  of  reaction.  Thus  the  higher  the  recy¬ 
cling  ratio  gets,  the  lower  the  electric  power  output  becomes. 
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system  design.  The  first  of  them,  the  comparison  of  differ¬ 
ent  reforming  concepts,  is  an  illustrative  pleading  in  favour 
of  the  DIR  concept,  the  direct  internal  reforming  in  the  an¬ 
ode  channel.  The  benefit  of  fuel  cell  cascades  compared  to 
single  cells  is  estimated  in  the  second  application,  proving 
its  usefulness  for  conceptual  system  design.  Results  indicate 
that  considerable  power  increase  can  be  expected,  but  the  ad¬ 
ditional  hardware  effort  might  spoil  the  benefit  for  smaller 
systems.  The  third  application  demonstrates  that  gas  recy¬ 
cling  spoils  the  system  performance.  This  application  also 
reveals  the  limitations  of  the  steady  state  anode  model,  as 
temperature  effects  are  not  considered.  Further  applications 
of  the  model  and  its  accompanying  diagram  are  possible. 


Fig.  13.  Trajectories  for  systems  without  and  with  anode  gas  recycling. 
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Appendix  A 

A.  1.  Derivation  of  the  steady  state  anode  model 


The  partial  and  the  total  mass  balance  for  the  anode  chan¬ 
nels  in  molar  terms  under  the  given  model  assumptions  read 
[7,8]: 

^“j7  -  Y  ~  X0>j)Dajrj 

j= ref, ox 

(A.l) 

Y  D/Da;G 

7= ref,  ox 

(A. 2) 

with 

>< 

II 

o 

II 

>< 

5’ 

(A.3) 
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cr 

ii 
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II 

X 

(A.4) 

The  steady  state  anode  model  is  a  simple  and  easily  un¬ 
derstandable  tool  for  the  analysis  of  various  aspects  of  the 
interaction  of  reforming  and  oxidation  reaction  in  high  tem¬ 
perature  fuel  cells.  It  is  small,  mathematically  simple,  fast 
to  solve,  flexible  and  can  easily  be  interpreted  by  physical 
means.  The  phase  diagram  of  the  model,  the  conversion  di¬ 
agram,  is  a  useful  and  intuitive  tool  for  the  display  and  in¬ 
terpretation  of  the  simulation  results.  Although  the  model  is 
applied  to  MCFC  here,  it  can  easily  be  adapted  to  fit  for  SOFC 
as  well.  The  basic  phenomena  shown  here  are  the  same  for 
both  fuel  cell  types. 

All  three  applications  which  are  shown  here  yield  qualita¬ 
tive  results,  one  of  them  gives  quantitative  results  suitable  for 


Expressing  the  molar  fraction  via  reaction  extents  (Eq.  (7)), 
these  equations  transform  to: 


d h 
d? 


Da/Q 

VHj,  max 


1  +  Y 

y'=ref,ox 


(A.5) 


d  y  _  Ej=ref,oxk/(d|;/dO£/,max  _  f'iK) 

df  1  +  X!  j=ref,ox  max  f{0 

with  the  boundary  conditions 


(A.6) 


m  =  0)  =  b  in  (A. 7) 

K(C  =  0)  =  rin  (A. 8) 
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The  solution  of  Eq.  (A. 6)  and  its  boundary  condition  in  Eq. 
(A. 8)  is 


y  =  r in  x 


1  +  X 

j= ref,ox 


(A. 9) 


Inserting  this  into  Eq.  (A. 5),  we  receive  the  desired  ODE 
describing  the  progress  of  the  extents  of  reaction  along  the 
channel: 


=  Da.;r;  x  1 

d?  Fn  £/,  max 


(A.  10) 


A.2.  Reaction  kinetics 


The  reaction  kinetics  is  described  using  potential  kinetics 
combined  with  Arrhenius  terms  for  the  temperature  influence 
and  Butler- Volmer  kinetics  for  the  electrochemical  reactions. 
The  molar  fractions  in  the  following  reaction  rate  expressions 
are  calculated  from  the  extents  of  reaction  according  to  Eq. 
(7).  The  reforming  reaction  rate  reads: 


r ref  =  exp 


X 


2  XCQ2Xh2  | 

Xch4Xh2o  Keq  iei(l})  j 


(A.ll) 


The  electrochemical  oxidation  rate  reads: 


rnx  =  exp  I  Arrox  I  -J-  -  ^  1  1  x  exp 


x  Xh2  -  Xh2oXco2  x  exp  -2 


A0a  A(/>oxo(^7) 
& 


,  A</>a— A0OX  0(ff) 
d 


(A.  12) 


The  electrochemical  reduction  rate  corresponds  to  the  de¬ 
scription  of  the  superoxide  mechanism  from  Prins-Jansen 

[9]: 


x  exp  I  2.5 


x  1  Xco 


'2,C 


x  exp 


red 

A  0c  —  A0red,o(^) 

-  y0-75  y"0-5 
a02,caco2,c 

,  A0C  A0red,o(^) 

i - 

$ 


(A. 13) 


Note  that  the  overpotentials  of  both  reactions  can  be  found 
in  the  argument  of  the  exponential  function,  as  is  shown  here 
for  the  anodic  reaction: 


(A.  14) 


Here,  the  Nernst  voltage  at  standard  gas  composition  occurs 
as  a  temperature  dependent  function,  A(/;ox,o(T). 

A. 3.  Cathode  gas  composition 

The  cathode  gas  composition  corresponds  to  the  exhaust 
gas  of  a  complete  combustion  of  the  feed  gas  with  a  specific 
amount  of  air.  The  air  amount  is  determined  via  an  air  number 
according  to: 

C-dir  =  - —  X  Tfeed2 XCH4 , feed  (A.  15) 

X02,air 

The  total  mass  balance  of  the  combustion  yields  the  cathode 
outlet  flow: 


rc  —  Xfeed  +  Xair 

(A.  16) 

The  partial  mass  balances  deliver  the  molar  fractions  of  the 
cathode  gas  as: 

r reed  XCH4 ,  feed 

XC02,c  = 

/  c 

(A.  17) 

^feed(XH-)0,feed  +  2xCH4,feed) 

XH20,c  =  r 

/  c 

(A.  18) 

^~airX02,air  -Oeed^XCH^feed 

X02,c  = 

/  c 

(A.  19) 

r airXN2,air 

XN2,c  = 

/  c 

(A. 20) 

There  are  no  further  components  present  in  the  cathode  gas. 


A.  4.  Used  model  parameter  values 
All  variables  are  dimensionless. 


Day  Arry  Equilibrium 

Reforming  25  84.4  ^eq,ref  =  exp  (26.22  — 

Oxidation  5  21.6  A0OX,O  =  28.26  -  19.84?? 

Reduction  0.3  31.2  A0red,o  =78.00  —  23. 06# 


Feed  Air/Others 
S/C  =  2.5  Aair  =  2.2 
/feed  =  1  XO2,air=0.21 

F=(l  +S/C)/8 
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